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Abstract 

We study bifurcations associated with stability of the lowest stationary point 
(SP) of a damped parametrically forced pendulum by varying (the natu- 
ral frequency of the pendulum) and A (the amplitude of the external driving 
force). As A is increased, the SP will restabilize after its instability, destabi- 
lize again, and so ad infinitum for any given loq. Its destabilizations (resta- 
bilizations) occur via alternating supercritical (subcritical) period-doubling 
bifurcations (PDB's) and pitchfork bifurcations, except the first destabiliza- 
tion at which a supercritical or subcritical bifurcation takes place depending 
on the value of ujq. For each case of the supercritical destabilizations, an infi- 
nite sequence of PDB's follows and leads to chaos. Consequently, an infinite 
series of period-doubling transitions to chaos appears with increasing A. The 

critical behaviors at the transition points are also discussed. 
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I. INTRODUCTION 



A damped parametrically forced pendulum (DPFP) with a vertically oscillating support 
is investigated. It can be described by a second-order non-autonomous ordinary differential 
equation (ODE) 

X + 2'K'~^x + 2tx{ijJq — A cos 27rt) sin 2t[x = 0, (1) 

where x is the angular position, 7 the damping coefficient, uo the undamped natural fre- 
quency of the unforced pendulum, and A the amplitude of the external driving force of 
period one. The overdot denotes the differentiation with respect to time, and all variables 
and parameters are expressed in dimensionless forms. 

The DPFP, albeit looking simple, shows a richness in its dynamical behavior. As the 
amplitude A is increased up to moderate values, transitions from periodic attractors to 
chaotic attractors and vice versa, coexistence of different attractors, transient chaos, and so 
on have been found numerically ||^-|^ and analytically 0. They have been also observed in 
real experiments [§,0. However, as A increases further, the DPFP exhibits new interesting 
dynamical behaviors not found in previous works, as will be seen below. 

Here we are interested in bifurcations associated with stability of the lowest stationary 
point (SP) with x = and i; = of the DPFP. The linear stability of the SP is determined 
by the linearized equation 

X + 2TT-fx + Att^{ujI - Acos2TTt)x = 0, (2) 

which is a damped Mathieu equation. For the undamped case with 7 = 0, stability properties 



of the Mathieu equation are given in |110| , pnj| . There exist an infinite number of disconnected 
instability regions of the SP in the ujq — A plane. These instability regions may be called 
"tongues", because their lower parts are of tongue shape (see Fig. 7-5 in Ref. They 
can be also labelled by an integer n, half of which correspond, in absence of forcing {A = 0), 
to parametric resonant values of Wq, i.e. = |. However, even a small amount of damping 
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leads to the presence of a non-zero threshold value At{n) of the amplitude necessary for the 
occurrence of the nth-order parametric resonnance [Q^- Moreover, At{n) grows rapidly 
with increasing n (see Fig. 100 in Ref. [0]). 

We first introduce the Poincare map for the DPFP in Sec. ^ and then discuss various 
bifurcations associated with stability of periodic orbits. With increasing A to sufficiently 
large values, the bifurcation behaviors associated with stability of the SP are investigated 
in Sec. |Tl| for a moderately damped case with 7 = 0.1. The damped Mathieu Eq. (§) has 
an infinity of alternating stable and unstable A ranges for any given uq- Hence, as A is 
increased, the SP undergoes a cascade of "resurrections" for any given c^o, i-e. it will resta- 
bilize after it loses its stability, destabilize again, and so ad infinitum. Its restabilizations 
occur through alternationg subcritical period-doubling bifurcations (PDB's) and pitchfork 
bifurcations (PFB's). On the other hand, the destabilizations occur through alternating 
supercritical PDB's and PFB's, except the first destabilization at which a supercritical 
or subcritical bifurcation takes place depending on the value of c^o- For each case of the 
supercritical destabilizations, an infinite sequence of PDB's leading to chaos follows. Con- 
sequently, an infinite series of period-doubling transitions to chaos appears with increasing 
A, which was not found in previous works. This is in contradistinction to the cases of the 
one-dimensional (ID) maps and other damped forced oscillators, for which only one single 
period-doubling transition to chaos occurs. In Sec. |I^, we study the critical scaling behav- 
iors at the transition points. It is found that they are the same as those for the ID maps. 
Finally, a summary is given in Sec. 0. 

II. STABILITY OF PERIODIC ORBITS, BIFURCATIONS AND LYAPUNOV 
EXPONENTS IN THE POINCARE MAP 

In this section, we first discuss stability of period orbits in the Poincare map of the 
DPFP, using the Floquet theory. Bifurcations associated with the stability and Lyapunov 
exponents are then discussed. 
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The second-order ODE ([l|) is reduced to two first-order ODE's: 



x = y, 



y = —2717?/ ~ ^tt^Uq — A cos 2nt) sin 2'7ix. 



(3a) 
(3b) 



The Poincare maps of an initial point Zo[= {x{0),y{0))] can be computed by samphng the 
points Zm at the discrete time t = m, where m = 1,2,3, ... . We call the transformation 
Zm Zm+1 the Poincare (time-1) map, and write Zm+i = P{zm)- 
The Poincare map P has the inversion symmetry such that 



SPS{z) = P for all z, 



(4) 



where z = {x,y), S is the inversion of z, i.e., S{z) = —z. If an orbit {zm} of P is invariant 
under S, then it is called a symmetric orbit. Otherwise, it is called an asymmetric orbit, 
and has its "conjugate" orbit S{zm}- 

We now study the stability of a periodic orbit with period q such that P'^{zo) = zq but 
P^{zq) 7^ Zf) for 1 < j < k — 1. Here P^ means the fc-times iterated map. The linear 
stability of the g-periodic orbit is determined from the linearized-map matrix DP'^{zq) of 
P'^ at an orbit point zq. Using the Floquet theory ||T^, the matrix DP'^ can be obtained by 
integrating the linearized differential equations for small perturbations as follows. 

Let z*(t) = z*(t+q) be a solution lying on the closed orbit corresponding to the g-periodic 
orbit. In order to determine the stability of the closed orbit, we consider an infinitesimal 
perturbation {6x{t),6y{t)) to the closed orbit. Linearizing Eq. (^ about the closed orbit, 
we obtain 



(5) 



where 
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An^ (uj^ - A cos 2nt) cos 2Tix*{t) -2n-f 



(6) 



Note that J is a 2 x 2 g-periodic matrix. Let W{t) = {w^{t),w'^{t)) be a fundamental solution 
matrix with W{0) = I. Here w^{t) and w'^{t) are two independent solutions expressed in 
column vector forms, and / is the 2x2 unit matrix. Then a general solution of the g-periodic 
system has the following form 



^ 6x{t) ^ 



5y{t) 



W{t) 



5y{Q) 



(7) 



Substitution of Eq. (|^) into Eq. (|^) leads to an initial-value problem to determine W{t) 



W{t) = J{t)W{t), W{0) = I. 



(8) 



It is clear from Eq. (J^) that W{q) is just the linearized-map matrix DP'^{zq). Hence the 
matrix DP'^ is calculated through integration of Eq. (|]) over the period q. 
The characteristic equation of the linearized-map matrix M(= DP'^) is 



- trM A + det M = 0, 



(9) 



where trM and detM denote the trace and determinant of M, respectively. The eigenvalues, 
Ai and A2, of M are called the Floquet stability multipliers. As shown in |13[, detM is 
calculated from a formula 



detM = e/o'f^'^*. 



Substituting the trace of M (i.e., trJ = —2n'y) into Eq. ([10|) , we obtain 



detM = e-^^^"?. 



(10) 



(11) 



Hence, the Poincare map P is a two-dimensional (2D) dissipative map with a constant 
Jacobian determinant (less than unity), like the Henon map ||14|| . 

The pair of stability multipliers of a periodic orbit lies either on the circle of radius 
£-1-79^ or on the real axis in the complex plane. The periodic orbit is stable only when both 
multipliers lie inside the unit circle. We first note that they never cross the unit circle, and 
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hence Hopf bifurcations do not occur. Consequently, it can lose its stability only when a 
multiplier decreases (increases) through —1 (1) on the real axis. 

A more convenient real quantity R, called the residue and defined by, 

1 + detM-t^ (12) 
2(l + detM) ' ^ ^ 



was introduced in [|15l to characterize stability of periodic orbits in 2D dissipative maps with 
constant Jacobian determinants. A periodic orbit is stable when < i? < 1; at both ends 
of i? = and 1, the stability multipliers A's are 1 and —1, respectively. When R decreases 
through (i.e., A increases through 1), the periodic orbit loses its stability via saddle- node 
or pitchfork or transcritical bifurcation. On the other hand, when R increases through 1 
(i.e., A decreases through —1), it becomes unstable via PDB, also referred to as a fiip or 
subharmonic bifurcation. For each case of the PFB's and PDB's, two types of supercritical 



and subcritical bifurcations occur. For more details on bifurcations, refer to Ref. [|T6| 



Lyapunov exponents of an orbit {z^} in the Poincare map P characterize the mean 
exponential rate of divergence of nearby orbits [|r^. There exist two Lyapunov exponents 
(Ji and (72 (cTi > (T2) such that (Xi + (72 = —2717, because the linearized Poincare map DP 
has a constant Jacobian determinant, detDP = e~'^'^"' . We choose an initial perturbation 
5zo to the initial orbit point zq and iterate the linearized map DP for 5z along the orbit to 
obtain the magnitute dm (= I^Zml) of ^Zm- Then, for almost all infinitesimally-small initial 
perturbations, we have the largest Lyapunov exponent ai given by 

ai = lim — In^. (13) 

m^oo m do 

If 0"! is positive, then the orbit is called a chaotic orbit; otherwise, it is a regular orbit. 

III. MULTIPLE PERIOD-DOUBLING TRANSITIONS TO CHAOS 

In this section, by varying two parameters uq and A, we study bifurcations associated 
with stability of the SP of the DPFP for a moderately damped case with 7 = 0.1. It is found 
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that with increasing A, the SP undergoes an infinite series of period-doubhng transitions to 
chaos for any given ujq. This is in contrast to ID maps and other damped forced oscillatiors 
with only single period-doubling transition to chaos. 

The stability diagram of the SP is given in Fig. 0. There exist an infinity of disconnected 
instability regions in the uq-A plane, which are separated by one connected stability region. 
The instability regions may be called "tongues" , because their lower parts are tongue-shaped. 
They can be also labelled by an integer n, half of which correspond to the parametric resonant 
values of uq (i.e., c^o = f ) in absence of forcing {A = 0) for the undamped case of 7 = 



|l|-y|,|TO|,|TT| . However, even a small amount of damping results in a non-zero minimal value 
At{n) of the amplitude necessary for the occurrence of the nth-order parametric resonance 
Furthermore, At{n) grows rapidly with increasing n (see Fig. p. Hereafter, each 
tongue of order n is denoted by T„. 

With increasing A, each tongue T„ is twisted to the left, and lies above the tongue T„_i. 
In such a way, tongues pile up successively, as shown in Fig. |1]. Consequently, there exist 
an infinity of alternating stable and unstable A ranges for any given ujo- Hence, as A is 
increased, the SP will restabilize after it loses its instability, destabilize again, and so ad 
infinitum for any given ujq. Such "resurrection" mechanisms are given below. 

Bifurcation behaviors at the tongue boundaries are investigated in details. They depend 
on whether the tongue-order n is odd or even. At the tongue boundaries of odd (even) 
order n, the residue of the SP is 1 (0). Consequently, PDB's and PFB's occur when tongue 
boundaries of odd and even order n are crossed, respectively. For example, the boundaries of 
Ti and T3 in Fig. |l] are PDB curves, while the boundary of T2 is a PFB curve. For the cases 
of PDB's and PFB's, there are two types of supercritical and subcritical bifurcations, which 
occur depending on where tongue boundaries are crossed. A saddle-node bifurcation (SNB) 
curve, at which a pair of stable and unstable orbits with period 2 (1) is born, touches each 
tongue boundary of odd (even) order at a boundary point [ujb{n), Af,{n)], and decomposes 
it into the supercritical and subcritical parts. As an example, see the three SNB curves, 
denoted by dash-dotted curves, touching the boundaries of Ti, T2, and T3, respectively. 



in Fig. |l|. On the lower left part of each tongue boundary, denoted by a solid curve, a 
supercritical bifurcation occurs. The remaining subcritical boundary curve starting from 
[ujb{n) , Ab{n)] first goes to the right, but it turns left at a point [utin), At{n)]. It consists of 
two types of subparts, denoted by short-dotted and dashed curves, on which a subcritical 
bifurcation takes place. On the subcritical segment with uJb{n) < ujq < ujt{n), the SP 
absorbes an unstable orbit born at a dash-dotted SNB curve and loses its stability. On the 
other hand, the stable orbit born by the same SNB undergoes an infinite series of PDB's 
leading to chaos. The accumulation points of such PDB's are denoted by open circles in 
Fig. |. 

When the SP loses its stability via supercritical PDB's and PFB's, the system is asymp- 
totically attracted to periodic attractors (born by the supercritical bifurcations) with the 
doubled period and the same period, respectively. However, for the subcritical bifurca- 
tion cases, the asymptotic states just after the instability of the SP may be periodic or 
chaotic, depending on which subparts of the subcritical boundaries are crossed. In Fig. |^, 
we fix different Uq and increase A to cross different subparts of a subcritical boundary of Ti. 
When a short-dotted boundary curve is crossed, the asymptotic state becomes periodic [see 
Fig. 2(a)], because the SP jumps to a periodic attractor born by an SNB after its instability. 
For this periodic case, with increasing A an infinite sequence of supercritical PDB's leading 
to small-scale chaos follows. However, when a dashed boundary is crossed, large-scale full 
chaos appears via intermittency [|1^, and hence the asymptotic state becomes chaotic [see 



Fig. 2(b)]. 

With increasing A to sufficiently large values, the bifurcation behaviors associated with 
stability of the SP are investigated in details for many values of uq- For a given c^o, the 
restabilizations of the SP occur via alternating subcritical PDB's and PFB's with increasing 
A, as shown in Fig. |I|. On the other hand, the destabilizations take place via alternating 
supercritical PDB's and PFB's, except the first destabilization at which a supercritical 
or subcritical bifurcation occurs depending on the value of Uq (e.g., for uq = 0.5 (0.65), 
the first destabilization occurs via supercritical (subcritical) PDB). For each case of the 
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supercritical destabilizations, an infinite sequence of supercritical PDB's leading to a pair 
of chaotic attractors follows and ends at a finite accumulation point. In each tongue, such 
accumulation points of PDB's, denoted by solid circles in Fig.l, seem to form a smooth 
critical line . Consequently, an infinite series of period-doubling transitions to chaos appears 
with increasing A. This is in contradistinction to the ID maps and other damped forced 
oscillators, in which only single period-doubling transition to chaos occurs. 

As an example of the multiple period-doubling transitions to chaos, consider the case 
ujq = 0.5. A bifurcation diagram along the vertical line uq = 0.5 is shown in Fig. |^. 
Through a supercritical PDB, the SP loses its stability at its first destabilization point 
A(i{l) = 0.100 218 ■ • •, and a symmetric orbit of period 2 is born. Unlike the case of the 
SP, the symmetric 2-periodic orbit becomes unstable by a symmetry-breaking supercritical 
PFB, which leads to the birth of a conjugate pair of asymmetric orbits with period 2. (For 
the sake of convenience, only one asymmetrical orbit of period 2 is shown in Fig. 3 [0.) 
However, as A is further increased, an infinite sequence of supercritical PDB's follows and 
ends at its accumulation point A^ (= 0.357 709 84 ■■■). The critical scaling behaviors of 
period doublings near the critical point Al are the same as those for the ID maps, as will 
be seen in Sec. 

After the period-doubling transition to chaos, a conjugate pair of small chaotic atttractors 
with positive largest Lyapunov exponent o"i appear. As A is increased, the different parts of 
a chaotic attractor coalesce and form larger pieces. For example, the chaotic attractor with 
(Ti ~ 0.091 shown in Fig. ^(a) seems to be composed of four distinct pieces for A = 0.3579. 
As shown in Fig. §(b), these pieces coalesce to form two large pieces with ai = 0.158 for 
A = 0.3582. However, beyond some critical point Ac{l) (— 0.3586), the chaotic attractor 
becomes unstable, and the system is asymptotically attracted to a rotational orbit of period 
1 born by an SNB. For A > ^4^,(1), the DPFP continues to exhibit rich dynamical behaviors. 
With increasing A, birth of new periodic attractors via SNB's, transitions from periodic 
attractors to chaotic attractors and vice versa, coexistence of different attractors, and so 
on are found until the SP restabilizes. (For more details on such dynamical behaviors, 
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refer to previous works However, with increasing A further, the DPFP exhibits new 

interesting dynamical behaviors not previously found. 

When the dashed subcritical boundary of Ti is crossed at the first restabilization point 
Ar{l) (= 3.150 509 ■■■), the SP restabilizes via subcritical PDB. An "inverse" process of 
the case of Fig. 2(b) occurs. There exists large-scale full chaos below ^^(1). When A 
increases through ^^.(1), the large chaotic attractor disappears, and the restabilization of 
the SP occurs with birth of a new unstable 2-periodic orbit. The residue R of the SP 
decreases monotonically from one, and becomes zero at the second destabilization point 
Ad{2) (= 3.224 230 ■ ■ •) on the supercritical PFB curve of T2. 

A second bifurcation diagram for uq = 0.5 is shown in Fig. ^ The SP becomes un- 
stable via symmetry-breaking supercritical PFB at its second destabilization point Ad{2), 
which results in the birth of a conjugate pair of asymmetric orbits with period 1. With 
further increase of A, a second infinite sequence of supercritical PDB's follows and ends 
at its accumulation point A2 (= 3.263 70315 ■ • •). The critical scaling behaviors of period 
doublings near A = A^ are the same as those near the first accumulation point A\. After 
the second period-doubling transition to chaos, a conjugate pair of small chaotic attractors 
also appears. They persist until some critical point Ac{2) (~ 3.263 862), beyond which the 
system is asymptotically attracted to an oscillating 2-periodic orbit born via SNB. As in 
the tongue of order 1, the DPFP exhibits diverse dynamical behaviors such as transitions 
between the periodic and chaotic attractors and the coexistence of different attractors in the 
region between Ac{2) and the second restabilization point Ar{2) (= 10.093 985 ■ ■ ■). 

When the dashed subcritical boundary of T2 is crossed at ^^(2), a subcritical PFB 
occurs. Consequently, the SP restabilizes with birth of a new unstable orbit of period 1. 
As A is further increased, the residue R of the SP monotonically increases, and becomes 
one at the third destabilization point ^4^(3) (= 10.097583 ■ ■ ■) on the supercritical PDB 
curve of T3. Since the order of T3 is odd, the subsequent bifurcation behaviors in T3 are the 
same as those for the case of Ti. That is, a third infinite sequence of supercritical PDB's, 
leading to a pair of small chaotic attractors, follows and ends at its accumulation point ^3 
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(= 10.099 660 93 ■ • ■). This third bifurcation diagram for uq = 0.5 is given in Fig. ^. The 
critical scahng behaviors of period doubhngs near are also the same as those near A^, as 
will be seen in the next section. 

We have also studied many other cases with different uo, and found multiple period- 
doubling transitions to chaos with increasing A. Such accumulation points are denoted by 
solid circles in Fig. |l|. In each tongue of order n, they form a smooth critical line A^{ujo)- 
Since the range of uq is < uq < u!b{n), each critical line of order n ends inside the tongue 
with order n. As mentioned above, a stable periodic orbit, born at a dash-dotted SNB curve, 
also undergoes an infinite sequence of supercritical PDB's. The accumulation points of such 
PDB's, denoted by open circles in Fig. |l], form another critical line. The two different critical 
lines joins at a point with uq = oj^in). Consequently, each critical line of order n extends to 
the outside of the tongue of order n. 

IV. CRITICAL BEHAVIORS OF PERIOD-DOUBLING BIFURCATIONS 

In this section, we study the critical behaviors (CB's) of PDB's for many values of ujq. 
The orbital scaling behavior and the power spectra of the periodic orbits born via PDB's as 
well as the parameter scaling behavior are particularly investigated. The CB's for all cases 
studied are found to be the same as those for the ID maps. 

As an example, we consider the case ujq = 0.5. The first three period-doubling transition 
points y4*'s [i = 1,2,3) are shown in Fig. |^. Only the CB's near Al are given below, 
because the CB's at the three transition points are the same. For this case, we follow the 
periodic orbits of period 2^ up to level k = 8. As explained above, for A = ^^^(l), the 
SP becomes unstable via supercritical PDB and a new symmetric 2-periodic orbit appears 
(see Fig. 3). However, the symmetric orbit of period 2 loses its stability by a supercritical 
symmetry-breaking PFB at A = 0.335 257 ■ ■ ■ . As a result, a conjugate pair of asymmetric 
2-periodic orbits appears. As A is further increased, each asymmetrical orbit with period 2 
undergoes an infinite sequence of supercritical PDB's, ending at its accumulation point Al. 
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Table | gives the A-values at which the supercritical PDB's take place; at Ak, the residue 
Rk of an asymmetric orbit of period 2^ is one. The sequence of converges asymptotically 
geometrically to its limit value A* with ratio 6: 

Aj. — A^^i 

= ^ ^ 5- (14) 

The sequence of 5^ is also listed in Table |. Note that its limit value 5 (~ 4.67) agrees well 
with that (= 4.669 • ■ ■) for a ID map Xm+i = f{xm) with a single quadratic maximum x* 
pO| . We also obtain the value of Al (= 0.357 709 845 3) by superconverging the sequence of 

For the ID map /, consider a 2*'-periodic orbit point x^'^^ nearest to the maximum point x* 
when the orbit becomes unstable. Then, the sequence of x^^^ also converges asymptotically 
geometrically to the maximum point x* with ratio a = —2.502 ■ ■ ■ ||2^. Note that the region 
near the maximum point x* is the most rarified region, because the distance between x*^^^ 
and its nearest orbit point f^'' ^(x'-'^^) is maximum. Hence, for the case of the Poincare map 
P, we first locate the most rarified region by choosing an orbit point z^^^ [= {x^''\ y^''^)] 
which has the largest distance from its nearest orbit point P^'' ^(2;^^^) for A = A^. The two 
sequences {x^^^} and {y'''^^} are listed in Table ||. Note that they converge asymptotically 
geometrically to their limit values x* and y* with the ID ratio a. 

_ y{k) _ y{k-l) 

"^''^ ^ ^{k+l) _ ^(fc) ^ ^y^^ ^ y{k+l) _ y{k) ^ (-'-^) 

The values of x* (= 0.091 126) and y* (= 0.735 292) are also obtained by superconverging 
the sequences of x'^^'^ and y'^^\ respectively. 

We also study the power spectra of the 2'^-periodic orbits {k = 1, . . . , 8) at the PDB 
points Ak- Consider the orbit of level k whose period is g = 2*^, {z^^ = (x^\y^^), m = 
0,l,...,g — 1}. Then, the jth Fourier component of this 2'^-periodic orbit is given by 



'iUJ 



1 m=0 



where cuj = 27ij /q, and j = 0, 1, . . . , g — 1. The power spectrum P^''\ujj) of level k defined 
by 
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P^'^\u,) = \z^^\u:,)\\ (17) 

has discrete peaks at = uoj. In the power spectrum of the next {k + 1) level, new peaks 
of the {k + l)th generation appear at odd harmonics of the fundamental frequency, uOj = 
27r(2j + l)/2('=+i) (j = 0, . . . , 2^= - 1). To classify the contributions of successive PDB's in 
the power spectrum of level k, we write 

p(^) = PooSiuj) + E ' E ' - (18) 

1=1 j=0 

where Pj^p is the height of the jth peak of the /th generation appearing at uj = ooij (= 
27r(2j + l)/2'). As an example, see the power spectrum P^^\uj) of level 8 shown in Fig. |7]. 
The average height of the peaks of the /th generation is given by 

<p''\^)=^'^i'pir- (19) 

It is of interest whether the sequence of the ratios of the successive average heights 

2/3W(/) =0(^')(/)/0('=)(/ + l), (20) 

converges. The ratios are listed in Table They seem to approach a limit value, 2/5 ~ 21, 
which agrees well with that (= 20.96 ■ ■ ■) for the ID map . 



V. SUMMARY 

Bifurcations associated with stability of the SP of the DPFP are investigated by varying 
two parameters c^o and A. As A is increased, the SP undergoes an infinite sequence of 
alternating restabilizations and destabilizations for any given uq. The restabilization and 
destabilization mechanisms are also given in details. A new finding is that an infinite series 
of period-doubling transitions to chaos appears with increasing A, which was not found in 
previous works. This is in contradistinction to the cases of the ID maps and other damped 
forced oscillators with only single period-doubling transition. The critical scalings at the 
transition points are also found to be the same as those of the ID maps. 
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TABLES 

TABLE I. Asymptotically geometric convergence of the parameter values for successive super- 
critical PDB's from an asymmetric 2-periodic orbit. 



k Ak Sk 

1 0.354163 288 011 

2 0.357 022 317174 5.286 

3 0.357563141135 4.692 

4 0.357678400 212 4.665 

5 0.357 703107 281 4.666 

6 0.357 708401983 4.668 

7 0.357 709 536 272 4.670 

8 0.357 709 779 136 



TABLE IL Asymptotically geometric convergence of the orbital sequences {xk} and {yk}- 



1 


0.094410 516 




0.719 956 679 




2 


0.088 901 931 


-1.933 


0.738 357 722 


-3.935 


3 


0.091750 680 


-3.085 


0.733 681 177 


-2.156 


4 


0.090 827396 


-2.261 


0.735 850 056 


-2.717 


5 


0.091235 660 


-2.635 


0.735 051829 


-2.398 


6 


0.091 080 705 


-2.436 


0.735 384 746 


-2.558 


7 


0.091 144315 


-2.538 


0.735 254 611 


-2.474 


8 


0.091119 256 




0.735 307 206 
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TABLE III. Sequence 2/3(*^)(/) [= + 1)] of the ratios of the successive average 
heights. 

k I 

3 4 5 6 7 

6 19.8 22.5 21.1 

7 19.8 22.1 21.2 21.5 

8 19.8 22.0 20.7 21.6 21.4 
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FIGURES 

FIG. 1. Stability diagram of the SP of the DPFP. There exist an infinity of "tongues" r„ of 

instability regions. For each tongue, a supercritical bifurcation occurs on the solid boundary curve, 
whereas a subcritical bifurcation takes place on the remaining dashed or short-dotted boundary 
curve. There are also SNB curves touching the tongue boundaries, which are denoted by the 
dash-dotted curves. The accumulation points of PDB's, denoted by solid and open circles, form 
critical lines. For other details, see the text. 

FIG. 2. Asymptotic states after the instability of the SP via subcritical bifurcations. A pair of 
symmetric stable and unstable orbit of period 2 are born via an SNB. The x-positions of the stable 
and unstable orbits are denoted by the solid and dashed curves, respectively. At a subcritical PDB 
point, the SP, whose a;-position is denoted by the dotted line, loses its stability by absorbing the 
unstable 2-periodic orbit. After its instability, (a) the SP jumps to the stable 2-perodic orbit for 
coo = 0.55, whereas (b) large-scale full chaos appears for uq = 0.6832. Note also that for each case, 
the stable 2-periodic orbit undergoes an infinite sequence of PDB's leading to small-scale chaos. 

FIG. 3. First bifurcation diagram for loq = 0.5. The SP2 and ASP2 denote the stable vl-ranges 
of the symmetric and asymmetric orbits of period 2, respectively. The PN also designates the stable 
A-range of the asymmetric periodic orbit with period N (N=4, 8, 16). 

FIG. 4. Chaotic attractors after the first period-doubling transition to chaos, (a) For 
A = 0.3579, the chaotic attractor with the largest Lyapunov exponent cxi ~ 0.091 is composed 
of four pieces, (b) These pieces merge to form two large pieces with ai ~ 0.158 for A = 0.3582. 

FIG. 5. Second bifurcation diagram for loq = 0.5. The ASPl and PN denote the stable 
A-ranges of the asymmetric orbit of period 1 and N (N=2,4, 8, 16), respectively. 

FIG. 6. Third bifurcation diagram for uq = 0.5. The SP2 and ASP2 denote the stable A-ranges 

of the symmetric and asymmetric orbits of period 2, respectively. The PN also designates the stable 
A-range of the asymmetric periodic orbit with period N (N=4, 8, 16). 
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FIG. 7. Power spectrum p(«) (w) of level 8 for A = As (= 0.357 709 779 136) 
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